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V ■ Abstract 

In this paper, a complete covariant quantization of generalized electrodynamics is shown through 
the path integral approach. To this goal, we first studied the hamiltonian structure of system fol- 
lowing Dirac's methodology and, then, we followed the Faddeev-Senjanovic procedure to obtain 
the transition amplitude. The complete propagators (Schwinger-Dyson-Fradkin equations) of the 
correct gauge fixation and the generalized Ward-Fradkin-Takahashi identities are also obtained. Af- 
terwards, an explicit calculation of one-loop approximation of all Green's functions and a discussion 
about the obtained results arc presented. 
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^ ■ 1 Introduction 

00 

The results that have been obtained for known theories using available theoretical tools are very 
| impressive: the agreement of QED^ with experiments, predictions of standard model and QCD4, and 
so many others. A point that warrants comment is the effectiveness of such theories up to a determined 
energy scale. Usually, a physics problem involves widely separated energy scales; this allows us to study 
the low-energy dynamics independently of the details of the high-energy interactions. The main idea 
is to identify those parameters that are very large (small) compared to the relevant energy scale of 
the physical system and let them go to infinity (zero). This provides a sensible approximation to 
the problem, which can always be improved by taking into account the corrections induced by the 
neglected energy scales as small perturbations. Effective theories constitute the appropriate theoretical 
tools to describe low-energy physics, where low is defined with respect to some energy scale. This idea 
of effective theories was proposed by Weinberg [TJ. 

The set of higher-order theories belongs to such effective theories. As it is known, the majority of 
physical systems described by Lagrangians depends, at most, on first-order derivatives. However, with 
the first development in formal aspects of higher-order derivative Lagrangians in classical mechanics 
by Ostrogradski [2], a new field of research was opened. 

The branch of higher-order derivative theories becomes very interesting, due to the fact that these 
additional terms are constructed in such way so as to preserve the original symmetries of problem. As 
a remark, it is important to say that this kind of theory has been shown to be a powerful method for 
consistent regularization of the ultraviolet divergences of gauge-invariant and supersymmetric theories 
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[3]. Also, the use of higher derivative terms becomes interesting regulatorr, by the fact that it improves 
the convergence of the Feynman diagrams [3J. 

More examples of systems treated with high-order Lagrangians that we can mention are: the 
study of the problem of color confinement on the infrared sector of QCD4 [5], the attempts to solve 
the problem of renormalization of the gravitational field [6], and a generalization of Utiyma's theory 
to second-order theories [7]. Although all these works improve the use of higher-order terms, the ones 
that most contributed to show the effectiveness of such terms in field theory was the contributions 
of Bopp jS], and Podolsky and Schwed [SJ, where they proposed a generalization of the Maxwell 
electromagnetic field. They wanted to get rid of the infinities of the theory, such as the electron 
self-energy (r _1 singularity) and the vacuum polarization current present on the Maxwell theory. The 
modification suggested by Podolsky and Schwed handle these unsolved problems and, also, gives a 
positive definite energy in the electrostatic case; also, as showed by Frenkel [TU], it gives the correct 
expression for the self- force of charged particles. In [7], it was shown that the Podolsky Lagrangian is 
the only possible generalization of Maxwell electrodynamics that preserves invariance under U(l). 

On theoretical and experimental framework, efforts have been made to determine an upper-bound 
value for the mass of the photon [TT] , the existence of a massive sector being a prediction of generalized 
electrodynamics. Along this line of thought, we believe that a way to set limits over Podolsky param- 
eter will be to study the Podolsky's photons interacting with standard model particles, and compare 
the obtained results with high-energy experiments. This idea and other purposes led Podolsky and 
some of his students to study the interaction of electrons with the Podolsky photons, which they called 
generalized quantum electrodynamics {GQED4) [12]. Among the points dealt with in their thesis, 
the most interesting was the calculation of electron self-energy at a one-loop approximation. They 
expected that the contribution of massive photons lead to a finite result; nevertheless, in the end, they 
found, as in the usual QED4, a divergent expression. Analyzing, now, the thesis results, we found 
a mistake in their treatment of theory, i.e., the choice of usual Lorenz gauge condition. However, 
this analysis was only possible due to the contribution of Galvao and Pimentel [13], which gives the 
first consistent quantization to Podolsky theory, where Dirac Hamiltonian formalism p3] was used 
with the correct choice of gauge condition, which they called the generalized Lorenz gauge condition. 
Also, they showed that, different from the usual Lorenz condition, the generalized one fulfills all the 
requirements for a good choice of gauge condition on the context of Podolsky theory. Indeed, one of 
the aims of this paper is to quantize GQED4, now, in the generalized Lorenz gauge condition. The 
Podolsky electrodynamics, by itself, takes account of several classical problems of Maxwell's theory, 
and it should be expected that the addition of Podolsky term into the QED4 Lagrangian with an 
appropriated gauge choice should give rise to interesting results. 

Based on all these facts pointed out above, we can conclude that higher-order theories deserve a 
deeper investigation. Therefore, this paper is intended to give a correct and transparent quantization 
of GQED4, the interaction of electrons with Podolsky's photons in four dimensional space-time. To 
improve our understanding of the features of the GQED4, we proceeded to calculate the radiative 
corrections of Green's functions. The main results of the paper will be closed formulas to the complete 
propagators and vertex function by functional methods (for a excellent review see [IS]), and it turns 
out that, with the correct gauge choice, the electron and vertex self-energy functions are finite at 
e 2 -order approximation. 

The work is organized as follow. In Sec. 2, we present a brief study of canonical structure of the 
theory and, then, construct the transition amplitude by the Faddeev-Senjanovic procedure [IB], which 



we believe is the most appropriated for our interest. In Sec. 3, we introduce the generating functional, 
which will generate all the Green's functions, photon and electron propagator, and vertex function; 
also, through it, we will derive the generalized Ward-Fradkin-Takahashi identities in Sec. 4. In Sec. 5, 
we evaluate and discuss the self-energy functions of theory at e 2 -order approximation. In order to 
avoid an awful reading, we place the most of calculation in the Appendices, and some useful identities, 
as well. Our remarks are given in Sec. 6. 



2 Transition Amplitude 

To construct the transition amplitude, we must first do a constraint analysis. Hence, before using the 
Faddeev-Senjanovic procedure, we will present a short study, showing the main points of Hamiltonian 
structure of GQED4: the evaluation of canonical momenta, followed by the determination of first- and 
second-class constraints, and, at last, the choice of an appropriated set of gauge conditions. However, 
it will be necessary to use the Faddeev-Popov-DeWitt method to get a convariant expression to the 
transition amplitude. Thus, we start with the Lagrangian density of GQED4, defined by0 



(1) 



C= l - Ud^ - - + e^AiP - jF^F^ + ^F^d? F a p, 

which, at classical level, is invariant under the local gauge transformations 

V/ (x) = e^H (x) , ^ (x) = A, (x) + U^X (x) . (2) 

In the Lagrangian (TjQ), we used the following definitions: the field-strengh tensor F Ufl = dyA^ — d^Ay 
and O = 7^0 M . The Lagrangian C preserves all symmetries of usual QED. The Euler-Lagrange 
equations following from the Hamiltonian principle with the corresponding boundary conditions are 

(id + eA-m\il> = 0, (l + a 2 a)d^F a ^ = eijjj a tp. (3) 

The canonical momenta, ix^ and <jr, conjugate to A a and T a , respectively, where T a = doA a are 
considered as independent variables, defined [T3], and given by 

* U735TTT " ^ k -^f— = F^ - a 2 W k d k d x F^ - flfa^H , (4) 



dc 2 



or 



V »°d x F 0X -OxF^ . (5) 



The canonical momenta associated with the fermion fields ip and ip are 

PA 55 dWM =! 2^^ (7) 

From the above momentum expressions, we shall study the constraint structure of the theory following 
the Dirac's approach to singular systems [2]. From equations dl])-© and the linear independence of 
constraints |14j . it is possible to obtain the following set of first-class constraints: 

Oi = <t> « 0, n 2 = 7r - d k( f> k wo, n 3 = d k n k + e^/ty « o, (8) 



1 We shall adopt, here, the metric convention rj^ — diag.(+, — , — , — ); the Greek and Latin indices runs from to 3 
and 1 to 3, respectively, and the spinorial indices are represented by capital Latin letters. 



XA = PA ~ 4 (7V) A -0, XA=PA-^ (^7°) A « 0, (9) 



and a set of second-class ones, 

where "~" represents the fact that the relations ([8]) and ([9]) are weak equations, according to Dirac's 
procedure. The constraint analysis presented here is justified by Faddeev-Senjanovic procedure to get 
the transition amplitude |16j . This point will become clear below. 

The transition amplitude in the Hamiltonian form is written in the following way 

Z 



D\i exp (% I d 4 x [vr^ (d Q A„) + <f> a (d T a ) - (doii) p - (d $) p-Hc\), (10) 
where the canonical hamiltonian Tic is given by 

Tic = 7r T° + TTjTi + + fa&To + 4>id k F lk - 3^ + ™H> - e^A^ + ^F kj F kj 

2 



+- (r, - d 3 A,Y - — [d^Y j - djffAoY . (11) 

The integration measure is defined by 

Dfi = D^DVD^DA^D^DpDpS (0,) det ||{fi a , det \\{ XA , Xb} b \\~ 1/2 (12) 

where = {f2, S,x, x}, and the functionals S are the gauge conditions that fix the first-class con- 
straints. Here, we will use the generalized radiation gauge condition 

Ei = r « 0, S 2 = A « 0, S 3 = (1 + a 2 U) 8 k A k » 0, (13) 



which as it is shown in [13], is an appropriated set of noncovariant gauge conditions for the first- 
class constraints (|8|). We notice that the determinant associated with the second-class constraints, 
det ||{xA;Xb}_bII) does not contain field variables, and so it can be absorbed in a normalization con- 
stant; we also show that the determinant between the first-class constraints ([8]) and the gauge fixing 
conditions (113p has the form 

det 1 1 £^11 = -(l + a 2 V 2 )V 2 . (14) 

Therefore, through the following manipulations-combining the equations (|12p and (|14p . substituting 
them into (|10p . and also carrying out momenta integrals and field variables- we find the following 
expression for the transition amplitude: 

Z = J DA^D4>Dij det \\-(l + a 2 V 2 ) V 2 \\5 (l + a 2 U) d k A k exp (i J d^xC^j. (15) 

Although equation f)15[) is correct, the noncovariant form is not good for calculation purposes. However, 
we can use the ansatz of Faddeev-Popov-DeWitt [T7] to achieve the desired covariant form for the 
transition amplitude. Then, choosing the generalized Lorenz gauge condition |13j 

n[A} = (l + a 2 a)d tl A ll = 0, (16) 

we finally obtain a expression for the covariant vacuum-vacuum transition amplitude 

DA^DipDip det j j - (1 + a 2 D) D| | x exp (i J d A x 

+ ^F^d a F^ - ± ((1 + a 2 D) <9^) z jj. (17) 

In this covariant gauge choice, we see that the Faddeev-Popov-DeWitt determinant does not con- 
tain field variables (the ghosts decouple from the gauge fields), and so, it can be absorbed into a 
normalization constant. 



i\) ( id - m + eA ) if) - -FavF^ 



3 Schwinger-Dyson-Fradkin Equations 



There are a lot of ways to extract the physical content of quantum field models, but the most ele- 
gant one is from the Green's functions using functional derivatives, which is a natural way to obtain 
such functions. The method of functional derivatives, which has been largely used by Schwinger, 
among others |181 119j . uses a generating functional from which all of Green's functions can be ob- 
tained by functional differentiation. These equations are also known as Schwinger-Dyson-Fradkin 
equations (SDFE), and the motivation to construct the SDFE's is the non-perturbative information 
that is provided for the theory. However, if we regard these equations only as a source of obtaining 
formal expansions in powers of the coupling constant, we shall obtain nothing new in comparison with 
perturbation theory. The problem of finding an effective method of solving those equations not based 
on perturbation theory is, at present, still far from any sort of satisfactory solution. 

In the present section, we will derive these relations for the photon and electron fields, and also 
for the vertex function. The first step is to define the generating functional 



2 bl,V,Jf J ]= / Dfi [tp, ip, Afj,) exp [iS e ff] 



(18) 



with the effective action given by 

S e ff = I (fix 



1 



where fj, r\ and are the sources [auxiliary mathematical device] for the fermion tp, anti-fermion tp 
and the gauge A^ fields, respectively. Let us stress that the components of fermionic fields [tJj, if) 
and their sources (rj, fj) are elements of the Grassmann algebra, and that A^ and its source J^, 
are c- numbers. From the generating functional (|18p . all the physical quantities of the theory can be 
obtained. Whenever possible, we will discuss the meaning of expressions of GQED4 and also its points 
of equivalence or inequivalence with the known results of QED4. 



3.1 Schwinger-Dyson-Fradkin equation for Photon Propagator 

We will derive and discuss, here, the properties of the complete expression of the gauge-field propagator 
in interaction with electrons. First, to obtain the corresponding photon SDFE we need to solve the 
following equation: 

SS eff 







SAp (x) 



S __6 5_ 

Sir) ' Sir) ' SiJn 



+ J M (x) 



(19) 



which, after evaluating the first term, can be written as 

1 



J» (x) 



□?7 



1 



'l + a z D 



5W 



+ re 



5W 



5W 



r 



5J V (x) Srj (x) 5rj (x) 



+ie- 



6 



5W 



5rj (x) \ 5rj (x) 



(20) 



The last equation represents the compact form of the non-perturbative equivalent to the Podolsky 
field equation, subject to an external source J M . The functional W present in (|20l) is the gener- 
ating functional for the connected Green's functions W [i], fj, J M ], which is defined by W [77, 77, J M ] = 



— i InZ [77, 77, Jn\. We also introduce the generating functional for one-particle irreducible (IP I) Green's 
functions V njj, ip, A^j through the Legendre transformation 



r [tp, i>, A^] = W [77, fj, J„] - / d 4 x ($7) + # + A^J^ 



(21) 



From the above definitions, we obtain expressions for (ip,ip, A^ in terms of (77,77, J^), and vice versa, 
being given by 



A -I— ti- - 1 6W l5W 

M i bJv ' i Sfj ' i Sr] 



Ju 



sr 



ST 



5T 



'SAP* 71 Sip' V Sip' 



Assuming the case that the fermionic sources are null, equation (j20|) is written as 



sr 



SAf, (x) 



+ 7(1 + ° 2n ) L ^ 



<l + a 2 U) UAa (x) + ie- 



5 



5W 



Sija (x) V fiv ( x ) 



where we have used the following set of projectors 



From identifying 



S 2 W[ V ,fj,J, 



□ ' 



y(x,y) = i 



(22) 
(23) 

(24) 

(25) 
(26) 



Srj (y) $V (x) 

as the complete electron propagator in an external field A^, which satisfies the following functional 
relation 



s 2 r 

i I d zS^bc (x, z) - = SbdS (x - y) 

owe (y) SipD (z) 



we can express (f2l"|) as 

sr 



SAf, (x) 



+ 1(1 + a 2 D) 



[1 + a 2 B) DAp (x) + eTr (7^ (x, x)) 



(27) 



(28) 



Now, differentiating (j28|) with respect to A u (y) and setting (x) = 0, yields 

5 2 r 



SA V (y) SAf, (x) 



+ 7(1 + « 2q ) ^ 
-ieTr 7^ 



□ (l + a 2 D) 5 (x - y) 

5 2 r N 



(29) 



SA V (y) \Sip (x)8ip (x) j 

The second term on right-hand side of (|29p can be evaluated immediately, giving a simple expression 



s 2 r 



SA v (y) \Sip (x) Sip (x) 



-1 



e I d ud w.y (u, x) r v (w, u; y) 5? (x, w) . 



(30) 



where we take into account the definition (|26p and have introduced the complete electron-photon 
vertex function 



er^ (x,y;z) 



s J r 



SAv (z) Sip (y) Sip (x) 



(31) 



A=V>=i/>=0 



Similar to the fermionic case ()27|) . the second derivative of V with respect to (x), gener- 

ates the inverse of the photon propagator <£>p V (x — y). From this fact, and substituting ([29]) into ([30]) . 
then follows the SDFE for the inverse of the complete photon propagator 



1 



V + j (1 + a 2 D) L 



/<// 



^l + a 2 D) nff(z-y) 



where the functional H^ u is known as photon self-energy function, and is defined as 
IV (x, y) = — ie 2 / d 4 ud 4 wTr (u, x) (x, w) T u (w, u; x) 



(32) 



(33) 



the (—1) factor comes from the fermionic loop in the usual way. The H^ u tensor describes the in- 
teraction of a photon with the electron-positron field, and this interaction consists of creation and 
annihilation of virtual pairs. Equation (|33|) in momentum representation assumes the form 



IV (k) 



ie 



(2tt)< 



d pTr 



y {p) 1^ (j> - k) T u (p, k; p-k) 



(34) 



From the expression ([32]) , we can compute the gauge- field propagator in a perturbative way, order- 
by-order, in the coupling constant e. An explicit calculation, analysis, and discussion at the lowest 
order of radiative correction of the Green's functions will be presented into Sec. 5. Indeed, since the 
bosonic IP I function and complete photon propagator satisfy the identity T ^ a {k)^ av (k) - 
is possible to find the following solution to the complete photon propagator: 



i%u (k) 



+ 



kuk-v 



k 2 [U(k) + (l-a 2 k 2 )} k 2 (l-a 2 k 2 ) 2 k 2 



(35) 



where II is called scalar polarization, which is introduced due Lorentz invariance of 11^, it has the 
structure 

IV (k) = {-Vnuk 2 + k^K) n (k) . (36) 

It should be noted that the expression (]35p shows that the n(/c) function is related with the transverse 
pole of photon propagator in momentum representation. The diagramatic representation of the SDFE 
for the photon propagator ([35]) is shown in the Fig. 1. 

The photon propagator at lowest order in perturbation theory, i.e., taking IT (/c) = on ()35|) . can 
be conveniently written as 



iD^u (k) 



r}^ ~ (1 



k 2 



1 



+ (1-20 



k^k v 



kfj,k u 



1 



k 2 - 

kfj, k v 



k 2 (k 2 



2 • 



(37) 



As it can be seen in (137p . the beauty of this expression is the appearance of the second term on right- 
hand side, which is originated from the Podolsky term and the generalized Lorenz condition. Note 
that this term has massive poles, m 2 = a~ 2 , which leads to a cancellation of IR divergences that are 
present in the first term, the Maxwell's term. Furthermore, the separation of massless (usually QED4) 
and massive modes in the propagator expression (]37|) in general gauge £ is owing to the linearity of 
fields in the gauge terms of the Lagrangian ([1]). By the relation between the Podolsky parameter and 
the mass of photons, it is possible to set a bound value for the photons mass, once we evaluate the 
parameter a [TTJ [20] . 



+ "WW 



n 



Figure 1: SDFE for the photon propagator 



3.2 Schwinger-Dyson-Fradkin equation for Fermionic Propagator 

In what follows in this subsection, we present the derivation of an integral expression to the complete 
electron propagator 5? . We also introduce the mass operator which contains all the radiative 
correction to the motion of electron (in the same sense of polarization operator II to the photons). We 
guide the derivation of SDFE for 5? in the same way as presented in last subsection for the photon 
propagator. We recall that the functional equation 



5S, 



eff 



_S 6 S_ 

6ir] ' Sifj ' SiJ, l 



r\ (x) 



5ip (x) 

•q (x) Z + i I id — m — ie^ 



Z[ti,ti,Jh] 
5 \ 5Z 



5Jv (x) J 5rj (x) 



(38) 



which is equivalent to the Dirac equation in presence of external sources, will define a relation between 
J? and jtft . Now, differentiating (|38|) with respect to r\ (y) and taking the fermionic sources going to 
zero, one gets 



i5 (x — y) = I id — m + eA(x) — ierf 



y(x,y) 



(39) 



SJt* (x) 

where we have used the definition (|26l) for 5? . Equation (139)) defines the non-perturbative connected 
two-point fermionic Green's functions. 

By means of a functional derivative identity, together with (|23p and ()30|) . and also taking the 
source going to zero, the last term of (|39p reads 



(40) 



y (x, y;A) = e / d 4 ud A zd 4 w% a (x - u) ,9> (x, w) T a (w, z; u) S" (z, y) . 



5J^ (x) 

Also, note that the electromagnetic potential A presented in the fourth term of (|39p vanishes in the 
absence of an external source; that is, A^ (x; = 0) = 0. Combining this fact with (I40p . the equation 
(1391) is rewritten as 



iS (x-y) = (i'fdp -m),y (x, y) - J d*zZ(x, z)& (z, y) (41) 
where the electron self-energy operator S introduced above is defined by the following relation: 

E(x -y) = ie V [ d 4 ud A w^ (x, w) % a (u - x) T a (w, y; u) . (42) 



(43) 



which, on momentum representation, is written as 

-,2 



(2tt)< 



^ / d^ky (p - k) % a (k) T a (p, k;p-k), 



If we denote conveniently YjS? (x,y; A) = J d 4 zY<(x, z)S? (z,y; A), then we can rewrite the equation 
ET|) in the following suitable form: 



(ij^d^ — m — £) y (x, y) = i8 (x — y) 



(44) 



Figure 2: SDFE to electron propagator 



Moreover, introducing the so-called mass operator 

^(x, y) = m5 (x — y) + y) (45) 

into (|44|) . we obtain that the complete electron propagator in momentum representation assumes the 
form 

- Jt (p) 7^p M -m-lu{p) 
which states the relation between the electron propagator and the mass operator. The SDFE corre- 
sponding to the electron propagator is presented in Fig. 2. Equations (|44p and f|46[) show that the elec- 
tron propagator is the Green's function for an equation similar to the Dirac equation (p — m — T,)^ = 0, 
but differing from the latter by the addiction to the bare mass m of the quantity S. For this reason, 
is called mass operator. 

In a similar way to operator II, we can say that the operator £ describes the interaction of the 
electron with its own electromagnetic field. This interaction consists in emission and absorption of 
virtual photons. 

3.3 Schwinger-Dyson-Fradkin equation for Vertex 

As it is already known [22], it is impossible to construct for QED^ a closed integral equation that 
expresses the vertex function V in terms of 5? and 3l and that, together with the equations (135p and 
(|46p . would give us a complete system of equations determining the Green's functions. Nevertheless, 
it is possible to find a relation connecting the vertex function T with 5? and ^ |21j ; however, different 
from other Green's functions, this relation contains only skeleton graphs [TS], i.e., connected graphs. 
But, for our purposes here, it is enough to consider this kind of approximation, due to the fact that, 
here, we have only interest in e 2 -order calculation. Thus, recalling the vertex function is formally 
obtained from 

5{y{x,y-A)r l 

eF(w) = 5A(z) (47) 



with 1 being the inverse of the fermionic propagator (|46p , the vertex function can be also decom- 
posed as 

T M 0, y; z) = -i'y^x - y)5{y - z) + A M (x, y; z), (48) 

where is denoted as the vertex part of graphs. The vertex function can be expressed in momentum 
space in terms of a new unknown quantity, the electron-positron kernel K, by means of an integral 
equation [21] 



v^p,p'-k) = -i lll 8{p + p' -k) + J ^L[iy( P > + q )r IM ( q + P ',p + q) 

X iy(p + q)]K(p + q,p' + q,q), (49) 

where p' and p are, respectively, the momenta of the emerging and incident electrons, while k = p—p' is 
the transferred momentum. K consists of graphs with two external electron and two external positron 
lines. Well, we have obtained, here, a closed integral equation for the vertex function; however, for 




Figure 3: SDFE for the vertex function 

practical calculations we did not accomplished much, because is expressed in terms of an unknown 
quantity - the kernel K. We shall write down the complete kernel K as a sum over skeleton graphs, 
which in first-order yields [21] 

iK(p,p',k) = (ie) 2 T»(p,p - k)% u (k)T u (p' - k,p'). (50) 

Therefore, from (|50p . we find out that the skeleton equation for the vertex function (J48]) written 
in Fourier representation is given by 

e 2 r 

r M (p,p'; k) = -i^6(p -p -k) + — -\I J d A qi<y (p' - q) (p - q,p - q; k) iS? (p - q) 

x r a (p ) p-q;q)i9 ua (q)T v (p',p' -k,k). (51) 

Figure 3 shows the vertex function. It is important to emphasize, here, that the operators S, II and 
A introduced above are functional of the Green's functions =5^, & and T, which means that the self- 
energy functions are coupled, and one the Green's function depends of another ones of lower order. 
Hence, we clearly see that this tower of equations is related. 



4 Ward-Fradkin-Takahashi Identities 

As it is well known, the generalized Ward-Fradkin-Takahashi (WFT) identities are, in general, identi- 
ties among Green's functions following from the existence of a symmetry. The goal of this section is to 
derive these gauge identities to GQED4. First, we will show the WFT identity satisfied by IP I gauge 
function, which leads to the transverse character of operator 11^. Next, we will derive the relation 
between the vertex function and the inverse of the complete electron propagator, which is known as 
the main WFT identity. At last, we will reproduce the main WFT identity in the k — > limit (null 
transferred momentum). The derivation of these identities is formally given as follows: starting from 
the generating functional (|18p and performing the infinitesimal transformations 



$ (x) =il)(x) + i\ (x) V (x) , 4 (x) = A» (x) + (x) , (52) 

and noticing that neither gauge fixing term nor the source terms are invariant under these transfor- 
mations, we find that the generating functional Z [rj, fj, J M ] satisfie the following equation of motion: 

i- (1 + a 2 U) 2 d^—^— - f?-f- + vj^r-, ~ \ <H 2 = 0. (53) 
L e£ Jfj, (x) orj (x) or] (x) e J 

The next step in deriving of WFT identities is to express (I53p in terms of the connected Green's 
functions W [i], fj, J M ] as 

- ~ f (1 + ^fc,™ - it™ + - V' = 0. (54) 

e£ ' J M (x) or] (x) or][x) e 



Finally, one can obtain the main quantum equation of motion for the theory by writing (|54p into an 
expression for the 1PI generating functional V N),ip, through (|2ip . Thus, one has the general 
equation 

- - (1 + a 2 U) 2 d»A»{x) - i-^-V ix) + (x) + -^TTTT = 0. (55) 

e£ dip (x) dip (x) e oA fJi {x) 

From equation (|55j), it is possible to derive all WFT identities. Thus, the first identity comes by 
applying the functional derivative of A u (y) in equation (f55l) at A^ = ift = ip = 

(x, y)-j(l + a 2 U) 2 d v 8 (x-y) = 0, (56) 

which, together with equation (|32p . imply in 

fc M n^ (k) = 0. (57) 

which shows the transverse character of operator IP 1 '. Now, the main gauge WFT identity follows by 
taking the derivatives of expression (|55p with respect to ip (y) and 4>(z) at A^ = tp = ijj = 0, which in 
momentum space takes the form 

k^ (p,p'; k=p-p')=i (2n) 4 [y- 1 {p) - (p - p')] , (58) 

with y- 1 the inverse of the complete electron propagator (|46p . 

Although the local gauge invariance at the classical level has been broken in the quantum theory 
through the gauge fixing procedure and source terms, the main WFT identity (|58p holds, inheriting 
its essence, without which the renormalizability cannot be guaranteed. 

In the limit of null transferred momentum, i.e., k — > 0, the equation (|58p leads to a relation 

^(p,p;0) = (2^ 4 A^-i( p ), (59) 

from this limit, also follows that 

A<>,p;0) = -/-£(p). (60) 



Both relations, (|59[) and (|60|) . hold in the same way that they do for QED4. 



5 Radiative Corrections of the Second Order 

In the preceding sections, we have derived integral equations to the Green's functions, the electron and 
photon propagators, and vertex function for the GQED^. Now, we will investigate the corrections 
to these functions in the first nonvanishing order of perturbation theory. The expression for the 
operator H^ u at e 2 -order does not differ from that of the QED4. This divergent result implies, in 
the same way as for the QED^ [22], the renormalization of electronic charge e and the introduction 
of renormalization constant Z3 in the GQED4. Although the electron self-energy function £ and 
vertex part A in e 2 -order are different from the usual corrections for the QED4, due to the presence 
of Podolsky's terms in the free photon propagator D^ u (|37p . the structure of divergences at this order, 
by power counting, remains the same as QED4, linearly and logarithmically divergent. At first glance, 
this fact seems to lead to infinity results for the other two self-energy functions of GQED4 at e 2 -order, 
thus, an explicit calculation of £ and A expressions becomes necessary. These calculations are also 
necessary to verify whether the main WFT identity (|58|) is still satisfied at this order. 



p — k 



Figure 4: Electron self-energy diagram 

In order to use the dimensional regularization procedure, the Lagrangian must have the right 
dimension (of internal loops); then, it is necessary introduce the t'Hooft mass fi. Thus, also considering 
the case £ = 1, we have 

C eff = ^d1>- - e//-^lV " ~F^ V + ^F^dvFl - - [(l + a 2 U) d^} 2 . (61) 

In the next two subsections, we shall compute the T, and A functions. We will show that both 
functions, £ and A, can be separated in two distinct contributions, the well-known contribution from 
the QED4 and a new one that we will call Podolsky contribution. However, we can observe by power 
counting that the Podolsky sector presents a divergent share with the QED sector; thus, we expect 
that they may cancel out the divergence of the GQED4. Now, we proceed to an explicit evaluation 
of the electron self-energy and the vertex part. 



5.1 Electron Self-Energy 

We begin by investigating the second-order electron self-energy function. This quantity corresponds 
to the diagram shown in Fig. 4. 

In accordance with equation (|43|) . the self-energy function £ can be written as 



where 



and 



E< 2 ) (p) = Yiqed (p) + S Pod (p) 

d d k 1 „ (p-k + m 



(62) 



'QED 



(p) = ¥ 



4-d 



(2ir) d k 2 



1 



(p — k) — m? 



^Pod (p) = ~ifJ> 



4-d 



d d k l A {P-k + mj y 



(2nf 



(p — k) — m? 



(* 2 ~ & 



k 2 k 2 1 



la 



k\ k v 



(63) 



The separation of YP^> in two contributions is only made possible by the linear structure of the 
free photon propagator (|37p . First, the regularized QED4 contribution for the electron self-energy is 
given by [22] : 

^QED (p) = g~2 ~ (P~ 4m) + Y.QED Finite (p) (64) 

with 

(1 — z) zp 2 — m 2 z 



1 1 f 1 

'QED Finite 0) = ~J^2 ^ ^ + l ) ~ 2m + X )] ~ ^2 J dz Wi 1 ~ Z ) ~ 2m \ ln 



where e = 4 — d is the dimensional-regularization parameter. 

Now, to evaluate the Podolsky contribution T,p oc i ([63]) . it is suitable to write it as 

£p 0d (p) = £4t(p) 

a=l 



47T/i 2 



(65) 



(66) 



so that the quantities E„ ' , are defined by 



J Pod 



(p) 



4-d 



rfrf/c T[p-k + m) 7a 



=Sd (p) 



A-d 



,(3) 
J Pod 



4-d 



(2vr) d 


(p — A;) 2 — m 2 
— + k 
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(2n) d 
d d k * 


(p — k) 2 — m 2 
(p — k + rnj k 


k 2 (k 2 ~ 


(2vr) d " 


(p — k) 2 — m 2 





(67) 



(68) 



(69) 



We are going now to calculate the expressions ^p od , (|6Tj) - (|69|) . To solve conveniently the momentum 
integration, we will use the Feynman parametrization and dimensional regularization. Using the both 
procedures in equation ([67]) . one can put it in the form 



10 J (27T) 



dz 



d d k ^ X [k-p-m)-fx 



where b 2 = (1 — z) 



zp 



(k - pz) 2 + b 2 

m 2 z. Introducing the change of variables k — > k — pz, we obtain 



4-d 



cfe7 [(l-z)p + m]7 A 



1 



/0 J (27T) d [k 2 + b 2 Y 

The k integration is carried out by using the identity (|1U2|) . so that (|70|) reads 

A-d / ^ 



J Pod 



(p) = ("I) 1 



(47T 



T r 2-- 

I V 2 



Gb 7 A [(1 - z)p + m] 7A [6 2 ]2- 



o 



Now, expanding ([7T]) around d = 4, we find that 



p ° dip) ^o"78^ (? - 4m) + S 



(1) M 

Pod Finite \P) 



(70) 
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where 



J Pod Finite 



(P) 



16vr 2 



[p(l+7) -2m(l + 2 7 )] 



l r 1 
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Jo 



(1 " *) (ZP 2 ~ 



m z 



(73) 



We can evaluate the other terms in a similar way; however, to avoid an extensive calculus, we present 
here only the results, leaving the explicit calculation of these quantities and other extensive expressions 
in Appendix B. The evaluated expressions of them are 



(2) (v) 

Pod Finite \P) 



Zpldip) = ~\^2 ( m -P)+^pld finite (P) 



(74) 
(75) 



with the finite parts given by ()113|) and ()117j) . respectively. 

Indeed, by combining the results of equations (|72p . (|74p and (|75p into equation (|66p . it follows that 
the regularized contribution of the Podolsky sector for the electron self-energy function is given by 



1 1 

Spod (p) = -g~2 - (P - 4m ) + ^Pod Finite (p) 



(76) 



where T,p od Finite is given by (|118j) . 

Therefore, it finally follows from a rearrangement of equations (|64p and ()76|) that the electron 
self-energy function (|62p . at e 2 -order, has the following expression: 



X (2) (P) 
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8vr 2 



16vr 2 
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16vr 2 7 jo 



d-z [(1 — z)p — 2m] In 
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dy [2m - (1 + 3y) p\ A 1 (p, x, y) 
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dy [(1 -y)p + m]p y A 2 (p, x, y) 



(77) 



with the quantities A\ and A 2 defined as 

(i - y)yp 2 - m 2 y - (l 



^i (p, y) = hi 
-4-2 (p, x, y) = 



(1 — y) yp 2 — m 2 y 
1 



1-y) yp 2 



m 2 y 



(i - y) (yp 2 - 



m^y 



Equation (|77p shows that the electron self-energy function E^, at e 2 -order, does not depends on /x, 
and that it is also free of divergences, which do not occurs in such ordinary QED^ as equation (I64p , 
This last feature is an interesting property of the theory. It seems that the Podolsky term in the 
Lagrangian f)61 1) acts like a natural regulator of the theory, due to its massive character. Nevertheless, 
a better analysis shows that the Podolsky term is not the only one responsible for the finiteness of 
electron self-energy in e 2 -order; the choice of the generalized Lorenz gauge condition (|16|) is also closely 
related to the finite result (I77p . Hence, we can conclude that the choice of the usual Lorenz condition 
to GQEDi leads to the divergent result for the self-energy of electron evaluated in the thesis advised 
by Podolsky p2]. 



5.2 Vertex Correction 

We now turn to the calculation of the vertex part A ll (p',p;q = p — p'^{p' \p) (|5T1) . where, as usual, p' 
and p are, respectively, the momenta of the emerging and incident electron, while q = p — p' is the 
momentum of the incident photon. The diagram that corresponds to this quantity is shown in Fig. 5. 

In the same way that occurs in equation (|62p for the electron self-energy function S, the vertex 
part ()50p also shows the splitting of its expression in two distinct contributions: 



A^) (p',p)=A" QED (p', P )+^ Pod (P^P) 
One contribution comes from QED4 



A QED (P'.P) 
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and another one from Podolsky sector, 
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Figure 5: Vertex part diagram 



The regularized QED^ contribution (|79p for the vertex part is known as 



A QED {P',P) = + MQED Finite {p'lP) i* 1 ) 



with A QED Finite (p'.P) S iven b y 

A" ( P > p) ~ — td X f 1 " du ^EMS 1 y 

^QED Finite [P>P)- 16n 2 J Q dX J Q d V A 2 ^ 



1 + -+ / / dyln 



o j o 



47T/i 2 

(82) 

where we have introduced the functions 

(p',p,x,y) = 6(1 - x - y)p , j fl p + 2mq u [j u , 7 M ] - 4(l-ac-y + 3xy)p.p'j IM + 2m 2 7 At 
+ 2a; (1 - a?) 7^p 2 + 2y (1 - y) (^) 2 7 M - 4y (1 - y)p (p'Y - 4x (1 - x)pp" 
- 4 (1 - x - y - xy) (pV + p(p') M ) (83) 

and 

A 2 = xp 2 (1 — x) + yp' 2 (1 — y) — Ixypp — m 2 (x + y) , (84) 

to simplify the notation of integrals. 

Again, as it has happened with the Podolsky contribution for the electron self-energy function 
£p 0( i, the vertex part Ap od f)80f) can also be written as three terms, 

a=l 

where we have defined each term in the following way: 

■ Ml) / r \ ■ 4-d [ d d k a ft -k + m p-k + m 1 

as (p', P ) - f ?-\ +m 7 m p-£+™ fc 1 (87) 

PodW f J (27r)d ( p /_ A;) 2_ m2 / (p-fc) 2 -m 2 fc2 ^2_^' ^ > 

A M3) / / \_ . 4- d [ d d k ~ ft -k + m p-k + m - 1 

7 (i^* - »)» - n^ 7 (p - W - m> W^W' 
To evaluate such integrals, we will proceed as we presented in the last subsection for the electron 
self-energy function £p 0( i- In this subsection, we will only calculate one term, equation (|86p . and 
present the results of another ones, (j87|) and (|88p . leaving the calculation of the last two terms in 



Appendix C. Also, we exhibit there some extensive expressions that appear throughout this subsection. 
Hence, recalling the Feynman parametrization, equation (j86|) can be expressed as 



Jo 



(89) 



(2n) d [p + A 2-^(l-x-y)] 3 
where we have replaced k by k — xp — yp' and defined the function 

(k,p' ,p,x,y) = 7 CT (1 — y) p — k — xp + m 7 M (1 — x)p — k — yp 1 + m j a 
for convenience. 

We now attend to the k integration of (f8"9"j) . Using the properties of the Dirac matrices (|l(J9p to 
separate the different k terms in the numerator and performing the momentum integration with the 
aid of equations (|102j) and (|103p . we find that 
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and 



(p ,p, x, y) = [(1 — y) p — xp — m] 7^ [(1 — x) p — yp — m] 
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(90) 



(91) 



and the measure 



1 rl-x 

d$ = I dx dy. 
Jo 



Equation ()90|) . expanded around d = 4, gives the following expression: 

A£J (P',P) = -~-f + Aj£> Fmtt e V,P) 

with 
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The evaluated expressions for Ap^j (p',p) and Ap^J are [see Appendix C] 



1 1 
1 1 



_y + a m(2) 

e ' ~ Pod Finite 



P',P) 



iv Pod yp >p) — 

Kfd(p',P) = - e ^ + A^ dFmUe ( P ',p) 
where the finite parts are given by (|123|) and (|127j) . respectively. 



(95) 

(96) 
(97) 



Therefore, when the equations (|94p . (|96p and (|97p are combined, we determine the regularized 
expression for the Podolsky contribution to the vertex part Ap od (|85D : 



(98) 



Substituting the results of the equations (|8ip and (|98p into the definition (|78p . we obtain that the 
vertex part A^ 1 at e 2 -order has the following expression: 



A Pod(p'>P) = -- € ^2^ + A Pod Finite (P'^P) 



where A Pod Finite is given by the equation (|128j) . 
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(99) 



where we have defined the functions Si,5]2 and S» as 
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(101) 



to simplify the notation of integrals, as we have been done with the functions 3 and A. 

As stated in the beginning of this section, we have shown that both radiative corrections, the 
electron self-energy and vertex part, are finite at e 2 -order. Equation (|99p indicates the independence 
of the vertex part with the t'Hooft mass /j,, and again, as it has happened with the electron self- 
energy function, the finiteness of is due to the Podolsky term plus the choice of the generalized 
Lorenz gauge condition. Another important point is that the finiteness of vertex part A^ and electron 
self-energy function £ implies that the main WFT identity (|58p is still satisfied at e 2 -order. 



6 Remarks and conclusions 

In this paper, the effects of the Podolsky term in the quantum theory of electron and photon interac- 
tions were analyzed. After a constraint analysis, the covariant transition amplitude was derived with 
the aid of the Faddeev-Popov-DeWitt ansatz in the generalized Lorenz gauge condition. The choice of 
this gauge is of great importance to the obtained results. Then, we proceeded by deriving the SDFE's 
of theory by functional methods, and three Green's functions have been determined: the photon Si^ v 



and electron =5^ propagators and the vertex function T^, equations ([35]) . (fl6|) and (f5T]) . respectively. 
Through these functions, we introduced the self-energy functions that contain the radiative corrections 
in all order in pertubation theory: the polarization tensor IL^, the mass operator and the vertex 
part A^, respectively. However, modifications of these expressions compared to the ones for QED4 
were observed only in mass operator and vertex part, resulting from the contributions of the Podolsky 
electrodynamic. Although such modifications are presented in all Green's functions, only the photon 
propagator (|3T|) presents changes at tree level. Moreover, the most interesting feature of this expres- 
sion is the fact that we could separate the usual contribution for the QED4 in a general gauge ^ from 
one that arises from the Podolsky theory, and that the IR divergences presented in QED4 terms are 
suppressed by the massive terms of Podolsky contribution. 

The derivation of WFT identities was also presented. The first identity (|56p showed that the 
transverse character of the polarization tensor Yi^ is also preserved in the GQED4 as in QED4. 
Immediately, we found the main WFT identity that relates the IP I vertex function and the complete 
electron propagator. The main WFT identity (|58p is responsible for holding the essence of the gauge 
symmetry in quantum level, without which the renormalizability of theory cannot be guaranteed. 

The last part of the article was devoted to the analysis of what the Podolsky contribution brings 
to the quantum theory at e 2 -order in perturbation theory. At this order of approximation, we verified 
that the photon self-energy function is divergent, showing that, if we claim to the renormalization 
theory, the eletronic charge needs to be renormalized. Now, for the other two corrections, interesting 
features appeared, and the free photon propagator performs an important role in these analysis, 
giving origin to a splitting of the correction functions in two distinct contributions: one from the 
usual QED4 and another from the Podolsky theory. This splitting makes it possible to study each 
contribution independently. Thus, since the QED4 contribution is well-known in literature, our task 
here was to calculate the Podolsky contribution to the electron self-energy function and to the vertex 
part. And, the obtained expressions for Podolsky contribution Sp 0( i and Ap orf , equations f76|) and 
([98]) . respectively, present the same divergent terms of the QED4, equations (f6l|) and (f8T|) . but with 
oppositive signs, showing, then, that at e 2 -order the X and A^ functions are finite. Although, here, 
we restrict ourselves to the case £ = 1 , these results can be generalized. It is possible to show that for 
£ 7^ 1, the divergences associated with the electron self-energy function and vertex part of the QED4 
are also canceled by the Podolsky contribution. And, as an immediate consequence of the finiteness 
of S and A^, we verified that the main WFT identity (|58p keeps being satisfied. 

As a final comment, the Podolsky parameter a, which appears in all the expressions evaluated 
here as a free parameter (as the inverse of photon mass), can have its range of values limited through 
applications of Podolsky theory. For example, we can evaluate now the physical quantity u(p')A^u(p) 
that is related to the form factors Fi(q 2 ) and i ? 2(? 2 ) of electric charge e, and to the anomalous 
magnetic moment of the electron, respectively. We expect to set a bound limit to the Podolsky 
parameter a through the use of precise experimental data from the electron magnetic moment, by 
calculating the form factor F2(q 2 ) for GQED4. This study is now under development. We can also 
express the quantum theory in a more formal and constructive method, through dispersion relations 
[23| . which can give more transparent results and, also, a direct evaluation of electron anomalous 
momentum. Another interesting issue is the study of the gauge properties of the propagators for 
the GQED4, constructing and analyzing the Landau-Khalatnikov-Fradkin transformation [21] for the 
theory. As mentioned before, a renormalization process for the photon propagator is necessary, due 
to the divergence present in the self-energy; although the divergence is the same as the QED4, the 



renormalization constant and, also, the running coupling constant may differ from the results for the 
QED^ due to the poles from the photon propagator expression (|37p . 

Going beyond of T = 0, we can study the GQED4 at finite temperature, and derive all the thermo- 
dynamical quantities of theory, including the energy-density distribution. And, following the idea of a 
recent study of the Podolsky electromagnetism at finite temperature [20], where a bound value was set 
to the Podolsky parameter a through the energy distribution using the cosmic microwave background 
radiation temperature, we can also use the cosmic microwave background radiation temperature to 
set a value to a through the thermodynamical quantities of the GQED4. These issues and others will 
be further elaborated ,the subject of deep investigations, and reported elsewhere. 
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A d Dimensions Identities 



As we made use of dimensional regularization in the evaluation of the radiative correction expression, 
we present here some useful d-dimensional identities associated with integrals, properties of gamma 
function, and Dirac matrices. 



A.l Integration in d-dimensions 

The useful results of integrals that appear throughout the paper are 

d d k 1 r («-f) 



/ 



(2yr) d (k 2 - m 2 ) a (4vr) | r(a)[ _ mT ^ 

d d k k^k v r(a-i-f) 



(2n) d (k 2 - m 2 ) a ~ 2 (4yr) f r (a) [_ m 2]«-i-f 
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A. 2 The Gamma Function 

An important property of the gamma function, with small e, is given by the following relation: 
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where 



ipx (n + l) = l + - + + --7 

Z n 



and 7 is the Euler-Mascheroni constant. We needed the formulae 



zT(z) = T{z + l), X~s ~ 1 - -In A? 
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as well. 



A. 3 Dirac Matrices 

The algebra of Dirac matrices in (i-dimensions is 



{7",7"} = 2»r 

where rf" v is the metric tensor in d-dimensional Minkowski space (with signature + 
Sfi = d; hence, 
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.), so that 
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In addition 

Tr (odd no. of 7 matrices) = 
Tr I = f (d) , Tr ^ a ^ a = f (d) r] c 

Tr 7^7 Vf? = / (d) [rf V 
where / (d) is an arbitrary well-behaved function, with / (4) = 4. 



B Calculus of E^ (p) and E®, (p) 

(2) (3) 

In order to evaluate the terms ^ Pod and ^ Pod , we will follow the same steps presented in the calculation 
of Sp] rf in the subsection 5.1. First, we recall the Feynman parametrization and the dimensional 
regularization. Thus, from (j68h . we obtain 
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(2vr) d [k 2 + 

where we have changed k k — py, introduced b 2 = (1 — y)yp 2 — m 2 + x-$, and used the equation 
(|93p for the measure cfc. Since the integral of the odd powers of k in numerator is zero, it is enough 
to evaluate the contribution of the even powers. Then, carrying out the k integration, the equation 
(|110P is written as 
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Indeed, expanding (lllip for d — > 4, we find that Ev , can be expressed as 
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The term Xp^ d (|69j) is evaluated following the same steps as in the previous ones through the Feynman 
parametrization and dimensional regularization, and also replacing k —> k — pz: 
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where b 2 = (1 — y) (yp 2 — 4y) —m 2 y, as we have defined in subsection 5.1. Carrying out the momentum 
integration now, we find for (|114p the expression 
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which in the limit cZ — 4 is written as 
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Therefore, from the results os equations (|73[) . (|113p and ()117|) we obtain the following expression 
of the finite part of the Podolsky contribution Sp d Finite- 
1 
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C Calculus of A^/ (?/, p) and A^/ (p' } p) 

In the same way as we did in Subsec. 5.2, we will proceed here into the calculation of terms Ap 2 j 
and Ap^j of the Podolsky contribution to the vertex part at e 2 -order. Now, recalling the Feynman 
parametrization and the dimensional regularization, the equation (j87|) is expressed as 
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where we have replaced k by k — xp — yp' and defined conveniently the function 
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After a manipulation of matrices 7 in (I119P and evaluating the momentum integration, one gets 
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where we have defined the functions Si, S2 and S3 and the measure d£ in (jlOOp and (|10ip . respec- 
tively. 



Now, the equation ()12ip in the limit d — > 4 assumes the form 
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with 
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Following the same steps as before, we find for A^S (p',p), equation (|88|) . the following expression: 



(123) 



<i d /c O m (k,p',p,x,y) 



(27r) d [ k 2 + A^-^(l-x-y)] 4 
Now, when we evaluate the momentum integration in (|124j) . we get 
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which in the limit cZ — 4 is expressed as 
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with the finite part written as 
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After a rearrangement of expressions equations (j95H . (|123|) and (I127D . we find that the expression of 



the finite part of Podolsky contribution A p 
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is written as follows: 
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